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Abstract 

We study the phase structure of the unpolarized and polarized two-flavor quark matters at zero 
and finite temperatures within the Nambu-Jona-Lasinio (NJL) model. We focus on the region, 
which includes the coexisting phase of quark-antiquark and diquark condensates. Generalizing the 
NJL model so as to describe the polarized quark matter, we compute the thermodynamic potential 
as a function of the quark chemical potential (//), the temperature (T), and the polarization 
parameter. The result heavily depends on the ratio Gd/Gs, where Gs is the quark-antiquark 
coupling constant and Gd is the diquark coupling constant. We find that, for small Gd/Gs, the 
"ferromagnetic" phase is energetically favored over the "paramagnetic" phase. On the other hand, 
for large Gd/Gs, there appears the window in the (//, T)-plane, in which the "paramagnetic" phase 
is favored. 

PACS numbers: 11.10. Wx, 11.30.Qc, 11.30.Rd. 
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I. INTRODUCTION 



Quantum chromodynamics (QCD) is the fundamental theory of strong interaction be- 
tween quarks and gluons. QCD is an asymptotically free theory and perturbation theory 
may be used at high density and/or temperature. One-gluon exchange interaction between 
two quarks in the color antitriplet channel is attractive. Then one expects that a quark 
matter undergoes a phase transition to color superconducting phases at high density and 
relatively low temperature. Around two decades ago, it was shown that such phase tran- 
sitions really take place at high density and low temperature Around a decade ago, it 
was found 2] that the color superconducting gaps are of O(lOOMeV), which are comparable 
with nuclear matter density. 

Quark matter is expected to be produced in heavy-ion-collision experiments, and might 
exist inside compact stars. If quark stars exist, they are quark matters in itself. Such quark 
matters are the system with intermediate density at low temperature, so that perturbative 
QCD is not applicable. Although lattice QCD is powerful for studying the system with 
vanishing net baryon density, it is not efficient, at the present moment, for studying nonzero 
baryon density systems due to the "sign problem". Then, for studying such systems, one 
should rely on some effective theories of QCD, such as the extended Nambu-Jona-Lasinio 
(NJL) model [3!]. The NJL model incorporates above-mentioned attractive nature of the one- 
gluon exchange interaction between two quarks in the color antitriplet channel. Furthermore, 
it successfully describes static properties of the pion (e.g., pion mass and decay constant) in 
the hadoronic phase through chiral phase transition (see, e.g., |4|). 

Theoretical studies using the extended NJL model as well as other approaches have 
disclosed possible existence of various color superconducting phases; regular 2SC phase, 
charge-neutral gapless 2SC phase, color-flavor-locked phase, and so on. Furthermore, it has 
been found that, at moderate baryon density region, the chiral and diquark condensates 
co-exist 0. (For recent reviews, see, e.g., @, Q].) 

Magnetic property of quark matter is yet another important issue. Tatsumi was the first 
who pointed out the possible instability of a quark matter against the spin polarization 
(ferromagnetism) 8|: Within the one-gluon-exchange approximation in QCD, he computed 
the energy density S of quark matter at zero temperature as a function of polarization 
parameter p, and found that, at relatively low density, £{p) decreases as p increases. Since 
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then, there appear quite a few papers that are devoted to possible ferromagnetism in quark 
matter [9j. 

In this paper, we generalize the Tatsumi's analysis [8| to the case of polarized two-flavor 
quark matters with moderate baryon density at zero and finite temperatures. We introduce 
two types of polarized quark matters; the one is the "magnetic-moment-polarized" quark 
matters and the other is the "spin-polarized" quark matters. For dealing with such systems, 
we generalize the extended NJL model and employ the mean-field or Hartree approximation. 

In Sec. II, we briefly describe the extended NJL model and generalize it so as to deal with 
polarized quark matters. In Sec. Ill, we derive the thermodynamic potential Q for the two 
types of polarized quark matters at finite temperature within the mean field approximation. 
(Concrete derivation is given in Appendix D.) We find that the thermodynamic potentials 
for two types of polarized quark matters are of the same form. Difference may arise when 
the residual interactions are taken into account. In Sec. IV, we present the results obtained 
through numerical analysis. As in various analyses within the extended NJL model, the 
results are sensitive to the quark-antiquark coupling constant Gs and the diquark coupling 
constant Go- We show that, in the case of Gd/Gs < 1.15, the polarized state is energetically 
favored. While, for Gd/Gs > 1.15, there appears the window in a (/i, T)-plane, in which 
the unpolarized state is favored, where fi is the quark chemical potential and T is the 
temperature. Sec. V is devoted to summary and conclusion. In Appendix A, we briefly 
review the role of the projectors (r = ±, s — ±), Eq. (12.5[) . In Appendix B, we show 
that the polarized quark numbers are conserved to the first order of coupling constants. In 
Appendix C, for completeness, we give the forms of the quark propagators in the polarized 
quark matters. 

II. EXTENDED NAMBU JONA-LASINIO MODEL AND ITS GENERALIZA- 
TION 

A. Extended Nambu Jona-Lasinio model 

For describing the two-flavor quark matters, we adopt the extended Nambu- Jona-Lasinio 
model with the scalar-, pseudoscalar-, and scalar diquark-channels taken into account, whose 
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lagrangian density reads |7| 

C = q(i$ - m ) q + G s [{qqf + {qil^rq) 2 ] + G D [(iq c ee b j 5 q) (iqee b 'j 5 q c )] , 

where mo is the current quark mass, q c = Cq T , q c = q T C with C = i r ) 2r y° the charge- 
conjugation matrix. The quark field is a doublet in a flavor space and a triplet in a color 
space, q = q ia with z = 1, 2 and a = r(ed), g(reen), b(lue). The Pauli matrices r = (r 1 , r 2 , r 3 ) 
and e = ir 2 act on the flavor space, while (e b ) a/3 = e a/% (e 7 * 5 * = 1) is a totally antisymmetric 
tensor in a color space. Although the coupling constants, G$ and Go, enjoy the relation 
Gd/Gs = 3/4, we regard, as in [7|, Gd/Gs as a free parameter. 



B. Generalization to the case of polarized quark matter 

For dealing with polarized quark matters, we first introduce a spin-polarization vector 
n M (p), which is obtained from its rest-frame form, n M (0) = (0, e z ) = (0,0,0, 1), through a 
Lorentz transformation 

= — , n(p) = e'+ ■ (2-1) 

m m (£ p + m ) 

Here E p (= a/p 2 + Wq) is the energy of the (anti)quark with momentum p. Obviously, we 
have the relations, 

n 2 (j9) = nlip) - n 2 (p) = -1 , E p n (p) = p ■ n(p) . (2.2) 
Using Eq. (12. ip . we introduce a set of projectors onto the spin-polarized states: 

V ± {p) = l -{l± 1 d{p)) ■ (2.3) 
We also introduce a set of energy projectors [3]: 

A± ® = W p [Ep± 70 mo)] • (2 ' 4) 

It should be noted that 7oA_(p) = (m /E p )A + (p) and 7 A + (— p) = —(mo/E p )A-{p), where 
A + (_)(p) is the projector onto the positive (negative) energy state. Now we introduce a set 
of projectors: 

pW(p) = V s (p)A+(p) = ~k+(p)V s {-p) 

pH(p) = V s (-p)A4p) = A4p)V s (p) , (2.5) 
P^P^p) = (T^PW® , (r, r' = ±, s, s' = ±) . (2.6) 



In Appendix A, we show that, when acting on the quark field q(x), P^(iV) projects out 
onto the "spin-up" ("spin-down") positive energy-state, while P±~^(iV) projects out onto 
the "spin-up" ("spin-down") negative energy-state. 
The projectors P^\p) enjoy the following relations: 

ToPf^To = P^hP), (2-7) 
75Pf } (p)75 = PL?(-£), (2.8) 
C(P^(p)) T C = -PW(p). (2.9) 

We are now in a position to describe the "polarized" quark matter. From the above 
observation, we introduce the "spin up" /"spin down" quark number, Q^± \ and the "spin 
up" / "spin down" antiquark number, , where 

Qi T) = J d 3 x Q« = J d 3 xq^(x)PP(iV)q{x) (r = ±, s = ±) . 
1. Magnetic-moment-polarized quark matter 

Since a "spin up" ( "spin down" ) quark and a "spin down" ( "spin up" ) antiquark feel the 
same electromagnetic force, we introduce two quark chemical potentials: fi + is conjugate to 
the net quark-number charge with "positive magnetic moment (MM)", 

Q + = Q<"> + Q w (= J d 3 x Q + } , (2.10) 

and /!_ is conjugate to the net quark-number charge with "negative MM" , 

Q_ = Q H + Q W ^= J £ x Q_^j . (2.H) 

It can readily be shown that [Q +) QJ\ = 0. We call the quark matter with /i + ^ the 
"MM-polarized" quark matter. For dealing with such quark matters, the term that depends 
on ji + and /i_ should be added to the Lagrangian density; 

I -> I + fi+Q+ + A«-Q- • (2.12) 

The polarized charges Q + and Q- are not conserved. We show in Appendix B that they 
are conserved up to and including the first order of coupling constants Gs and Gd- It is to be 
noted in passing that, when the mean-field approximation is employed, even the unpolarized 
charge Q + + Q_ turns out not to be conserved 1 . 



1 We thank A. Oguri for discussion on this point. 
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2. Spin-polarized quark matter 



We define "spin-polarized" quark matters by introducing two chemical potentials, //+ and 
/i + is conjugate to + Q^" 1 (= Q' + ) and /i„ is conjugate to Q_ + (= Q'J)- 
One can show that [<3V?Q-] = 0- m this case, the counterpart to Eq. (I2.12p is 

£ ^ £ + fi+Q! + + fi-&- . (2.13) 

We see in Appendix B that the charges Q' + and Q'_ are conserved up to and including 
0(G s ,G D ). 

C. Mean field approximation 

We employ the mean field approximation to get 

I = q{i$ -m)q- - A* b (iq°ee b j 5 q) 
a 2 A* b A b 

where {qi^^rq) = has been assumed p| and 

m = m + cr = m — 2Gs(4>4>) > (2-15) 
A" = -2G D (zq c ee b l5 q) , A* b = -2G D (iqee b l5 q c ) . (2.16) 

m in Eq. (12.15P is the constituent quark mass, which has received the contribution a from 
the chiral condensate, if any. Then, when we use in the sequel various quantities defined in 
the last subsection, we replace the current quark mass mo with m. 

From Eq. ( 12. 1411 . we see that the red and green quarks participate in the diquark con- 
densate, while the blue quarks do not. 



l -A b (tqee b l5 q c ) 



(2.14) 
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III. THERMODYNAMIC POTENTIAL IN THE MEAN-FIELD APPROXIMA- 
TION 

A. Magnetic-moment-polarized quark matter 

Throughout in the sequel of this paper, we use the imaginary-time formalism. The grand 
partition function is defined by 

Z = j P^Pgexp jjf 9 dr j d 3 x ^ + J^ s oj j , (3.1) 

where j3 = 1/T is the inverse temperature, r = ix° (0 < r < 0) is the Euclidean time, and 
C E = C\ 

^ I XO = —IT 

Since the exponent of Eq. (13. ip is bilinear in fields, Z is factorized as 

Z = Z coas x,Z b Z r g . (3-2) 

^const reads 

b ^ t = -^(— + — j, (3.3) 
where is the volume of the system. The contribution from the blue quarks is 

Z h = J Vq b Vq b exj>y^ dr J d 3 x q b + ±qf (Gq ) _1 </) } , (3.4) 

where [G^d)]- 1 = ~ C ([Gg (-i9)] _1 ) T C and then 

(Co ) 1 = - m ± to /^T-Pi T) (^V) . (3.5) 

For the contribution from the red and green quarks, Q = q rg , we have 



Z r g = / VQVQexp{ / dr / d 6 x 



3, 







+\Q C (GoT 1 Q c + \Q^-Q C + ^Q C A + Q 



(3.6) 



where 

A~ = -iAee fe 7 5 , A+ = -iA*ee b -f 5 . (3.7) 
We employ the Nambu-Gorkov formalism through introducing 



Q 



c 



Going into a momentum space, we have (cf. Eq. (13.51) ) 

( G ^ = E ^° + ^ ± ^ 7 Pi T) (-p1 , (3i 



T, S=± 

where 



■7 



Z7T 

Po = ip OE = j(2n + l) (n = ...,-2,-1,0,1,2,...). (3.9) 
Eqs. (13 .41) and (13.61) turn out, in respective order, to 

Z fe = / V^ b exp | ~ ^ (/3G x ) M/ b 1 = [Det (-(3G, 1 )} 1/2 , 

Here 

=((<*>"' V G -=f( G « _l1 ^ ^ |, (3,0) 

V o (Go-)"!// V A+ K) M; 

where (l/)*- 7 = 5 i - 7 in the flavor space, and (l^) a/3 = 5 a/3 — 5 ab 6^ b in the color space. 

Although the forms of the propagators, Go and G, are not necessary for our purpose, we 
give them in Appendix C for completeness. 

Zb is obtained from Z T ^ g by taking the limit, A = A* — > 0; 

. (3.11) 

Z A = A* = 

Then we concentrate on evaluating Z rj9 ; 



h\Z b = -hiZ r n 
2 



lnZ r , 9 = ^lnDet(-/?G- 1 ) , (3.12) 



the computation of which is given in Appendix D. 



Thermodynamic potential 



From Eqs. (JH2D, flESD, flEID (cf. Eq. flTX3|) ). and Eq. fl3TTB . we obtain for the thermo- 
dynamic potential, 



T 

O = -— InZ 

^ 2 lAI 2 



(7 



+ 



4G S 4G 



Pi f=± 



+ i{^) + 2Tln(l + e-^)}^ o 



(3.13) 
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where Nf = 2 is the number of flavor- degrees of freedom, and 



E v + a 



yU+ + /i- 



IAI 



1/2 



+ P 



(3.14) 



The three momentum integral in Eq. (13.131) is ultraviolet divergent, and, as usual, we 
introduce a momentum cut-off parameter A. The zero-temperature limit (T — > 0) of Q 
reads 

a 2 |A| 2 Ar [ dh 

P, cr=± 



n(T = o) 



4G S 4G S 

in Eqs. (13 . 13[) and (I3.15P yields the gap equations: 

(A = |A| 



(3.15) 



on on 

0, — r~~ — o 



(3.16) 



dm OA 

A solution (m, A) to Eq. (13.161) is a local minimum point of Q, when the following conditions 
are met: 

d 2 Q d 2 Q ( d 2 Q \ 2 



dm 2 d\A 2 \ \dmdA ) 
d 2 Q d 2 Q . 
9m 2 <9|A| 2 



> 



(3.17) 



Polarized quark-number density 

Let n^' be the "positive (negative) magnetic moment" net quark-number density, n^' 
is computed from Q through partial differentiation with respect to /i±: 

dtt 



n (±) (/i+,/i_,T) = 

d/i± 

n^ + ,^-,T) = N f J 70^(A-P)E 



O'Ep a 



Pi cr 



^ P % + |A| 2 



± p (fig) 



(3.18) 



6»(A - p) [n F (E p - /i ± ) - n F (fi p + (jl±)] 



(3.19) 
(3.20) 
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where E v ^ a = E p + aji (p, = (//+ + /x_)/2) and Tip 0*0 = \/{eP x + 1) is the Fermi distribution 
function. In the limit T — > 0, turns out to 

?3„ 



ri?}(T = 0)=N f J±L0{A-p) 



+ 



fi- E p 



-0 



y/(ji + E p )* + \A\* ^-E^ + IAP 



9\Ep-fL + 



AP 



x 



fi- E p 



where 5/x = /i + — //_. 



v/(/i-E p ) 2 + |A|2 



Tl £ 



|A|2-£„ 



(2tt) 



-e(A- P )9(fi ± -E, 



V) i 



(3.21) 
(3.22) 



B. Spin-polarized quark matter 

The counterpart to Eq. (13.81) is 

(^r 1 = E (Po + rE P ± 7oPi r) (-p) • (3-23) 

T, S=± 

Q is computed in a similar manner as in Appendix D and we obtain the same form as above, 
i.e., Eq. (ETT5I) with Eq. (1331) . 



IV. NUMERICAL ANALYSIS 



In this section, through numerical analyses, we will analyze the phase structure of the 
unpolaried and polarized quark matters at zero and finite temperatures. How the chiral 
condensate, diquark condensate, and the degree of polarization compete with one another 
will be discussed. 

We first fix the parameters of the model as in [zj]: The current quark mass m = 5.5 
MeV, the three- momentum cut-off A = 0.637 GeV, and the coupling constant Gs = 5.32 
GeV -2 , which reproduces the properties of pion in vacuum. As mentioned above, although 
the standard value of Gp/Gs is 3/4, we leave it as a free parameter. 
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A. Unpolarized quark matter 
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Fig. 1: Plots of m and A as functions of the chemical potential \x for T = 5fi = 0. 
(a), (b), (c), and (d) correspond to Gd/G s = 3/4, 1.0, 1.2, and 1.5, respectively. 



Throughout in the sequel of this section, values of all dimensionful quantities are given in 
unit of MeV. In jl], the two gaps m and A(= |A|) at T — are displayed for various values 



of Gd/Gs- For completeness, we reproduce them in Fig. 1. The panel (a 
the standard value Gd/Gs = 3/4, while (b), (c) and (d) correspond, as in 



corresponds to 
J, to Gd/Gs = 



1.0, 1.2, and 1.5, respectively. Close observation on Fig. 1 is given in [7], which we briefly 
recapitulate here. 

Fig. 1 (a) {Gd/Gs = 3/4) shows that the chiral phase transition and the color su- 
perconducting (diquark) phase transition take place at nearly the same chemical potential 
fi x ~ /ia — 330. These two phase transitions are of first order. There remains a small chiral 
condensate in the color superconducting (CSC) phase /i > which is a relic of the explicit 
chiral symmtery breaking (mo 7^ 0). 
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We define //a, as usual, to be the value of the chemical potential at which the diquark 
condensate starts to appear. We define /i x to be the chemical potential at which m starts 
to decrease from nearly the constant value (rather than the chemical potential at which the 
chiral symmtery is restored). From Fig. 1 (b) - (d), we see that /ia — A* x and, as Gp/Gs 
inceases, /ia — fi x decreases. The chiral condensate in the coexistence phase is due to the 
dynamical symmtery breaking. 

Here we like to add one finding. Within the accuracy of numerical computation, the value 
of ji at which the quark-number density, n, vanishes coincides with /iA, i-e., n(/z,T = 0) = 
for n < //a- 

2. TV 

For the shake of copmarison to the case of 5fi ^ in the next subsection, we display here 
the results for T ^ and 5fx = 0. 




100 200 300 400 500 600 100 200 300 400 500 

\i (MeV) \i (MeV) 



Fig. 2: The same as Fig. 1 for T = 50. 

In Fig. 2, we plot m and A at T = 50 for the same values of G D /G S as in Fig. 1. Fig. 
2 (a) shows that the value of /ia is nearly the same as in the T = case, and ji x < //a- 
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Both phase transitions are of second order. A in the region /x > //a is smaller than that 
in the T = case. Comparing Fig. 1 (b) - (d) with Fig. 2 (b) - (d), we can make similar 
observation, besides that A's in the region \i > fj,& at T = 50 are not appreciably smaller 
than their T = counterparts. 

In the following we refer the region where A = (m = O(m )) to as the hadronic-phase 
(CSC-phase) region, while the region in between them simply to as the double-broken phase. 
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Fig. 3: m, A against T for (a) {G D /G s ,n) = (3/4,350) and (b) (G D /G s ,n) = (1-5,200). 

Let us see the T dependence of m and A. In Fig. 3, choosing the double-broken phase 
region, we plot m and A against T. The panel (a) [(b)] corresponds to Gjj/Gs = 3/4 [1.5] 
and /z = 350 [200]. We see that A decreases as T increases and becomes zero at the critical 
temperature Ta. Above Ta, the color normal phase is realized. With respect to the chiral 
symmetry, the panel (a) shows that, as T increases, m starts to decrease at T ~ Ta, and the 
chiral symmetry is restored gradually. On the other hand, in the case of panel (b), m starts 
to decrease from the beginning (T = 0) and the chiral symmetry is restored gradually. We 
have observed the similar behavious for other values of fi. 

Finally, in Fig. 4, we display the phase diagrams. The panels (a), (b), and (c) correpond 
to Gd/Gs = 3/4, 1.2, and 1.5, respectively. 

B. Polarized quark matter 



We have, among others, three parameters fi + , and T. We difine \x through 



H+ = M- 



- [n (+) (/i + ,/i_,T) +n ( ) (/i+,/i_,T 



(4.1) 



where n^(/x+, T) is as in Eqs. (13.181) - (I3.20p . The procedure of numerical computation 
for various quantities goes as follows. 
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Fig. 4: Phase diagrams for G D /G S = 3/4 (a), 1.2 (b), and 1.5 (c). 
H, CSC, and DB stand for hadronic, color-super-conducting, and double-broken phases, 

respectively. 



1) We first fix T, fi(> 0), and //_(> 0). 

2) Then, we solve two gap equatios, Eq. (13.161) . and Eq. (14.1 p simultaneously, and obtain 

the solution for m, A), and compute £/i = /i + — 

3) Finally, we check if the two conditions (13.171) are met. 



1. 5u dependence of m and A 

Double-broken phase: We first study the parameter region where the unpolarized quark 
matter is in the double-broken phase. In Fig. 5, we plot m and A agaisnt 5fi for different 
values of (Gd/Gs, T, fi): The panels (a), (b), (c), and (d) correspond to (3/4,0,400), 
(3/4,50,400), (1.5,0,300), and (1.5,50,300), respectively. Fig. 5 (a) shows that, as 6fi 
increases, the first-order phase transition occurs at 5fi = 5fi p ~ 220. For <5/i > 5fi p , the 
system is in the color normal phase. At 5fi = 5fi p , a gap appears in m. Fig. 5 (b) shows 
similar behaviors as Fig. 5 (a) , but no gap appears in m at 5fi = 5fi p . Figs. 5 (c) and (d) 
shows that, in the case of Gd/Gs = 1-5, both m and A are nearly independent of 5fi, and no 
phase transition occurs. We have carried out the same computation for different parameter 
values and found qualitatively the same results. 

Hadronic phase: We study the parameter region where the unpolarized quark matter is 
in the hadronic phase. In Fig. 6, m and A are plotted against Sfi for (Gd/Gs, T, /i) = 
(1.0, 50, 100) (cf. Fig. 2 (b)). We see that m and A (= 0) are nearly independent of 5/i and 
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Fig. 5: m and A against S/j,. (a), (b), (c), and (d) correspond to (Gd/Gs, T, /x)=(3/4, 0, 400), 
(3/4, 50, 400), (1.5, 0, 300), and (1.5, 50, 300), respectively. 
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Fig. 6: m and A against 5/j, for (Gd/Gs, T, Fig. 7: m and A against <5/x for (Gd/Gs, T, 



=(1.0, 50, 100). 



=(1.0, 50, 600). 



the quark matters are in the color normal phase for all 5fi. We have carried out the same 
computation choosing different parameter values and found qualitatively the same results. 
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CSC phase: Finally, we study the parameter region where the unpolarized quark matter is 
in the CSC phase. In Fig. 7, we plot m and A against 5/j, for (Gd/Gs, T, fi) = (1.0,50,600) 
(cf. Fig. 2 (b)). As 5ji increases, a phase transition occurs at 5/j, = 5/i p and the quark 
matter turns out to the color normal phase, while (small) m is nearly independent of 5/j,. 
We have carried out the same computation choosing different parameter values and found 
qualitatively the same results. 

2. Thermodynamic potential 

Here we compute the thermodynamic potential Q(/j,T; 5/j) as a function of 5/i, for dif- 
ferent values of (//, T), and discuss their physical implications. 

Two typical behaviors of Q are shown in Fig. 8. In the system with Q as in the panel 
(a), the "polarized" or "ferromagnetic" phase is realized, and we refer such an Q to as the 
"ferromagnetic (FM)" type. On the other hand, in the system with f2 as in the panel (b), 
the unpolarized phase is realized, and we refer such an f2 to as the "normal (N)" type. 




5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 40 

8n (MeV) 8n (MeV) 

Fig. 8: fi(T=0) against 5/i for /j = 400. The panels (a) and (b) correspond to Gd/Gs = 3/4 and 
1.2, respectively. 

Through numerical computation of Q, we have found the following results. 

1) For Gd/Gs < 1.15, f2 is of the FM-type in whole region of (/i, T). As an illustraion, we 

depict in Fig. 8 (a) = 400, T = 0; 5p) as a function of 5p, for Gd/Gs = 3/4. 

2) For G D /G S > 1.15, there appears a window in the (/i, T) -plane, in which f2 is of the 

N-type. We show in Fig. 9 the phase diagrams in the (/j, T) -plane. The panel (a) 
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Fig. 9: Phase diagrams with respect to polarization. The panels (a) and (b) correspond to 
Gd/Gs = 1-2 and 1.5, respectively. 

FM and N stand for the "ferromagnetic" and "normal" phase, respectively. 



[(b)] corresponds to G D /G S = 1-2 [1.5]. As an illustration, we plot in Fig. 8 (b) 
= 400, T = 0; 5/j,) against 5/j, for G D /G S = 1.2, which is of the N-type. 

To see how changes through the transition region from the FM-type region to the N-type 
region, we display Q{fi = 330, T = 0; 5fi) and n(jj, = 340, T = 0; Sfi) for G D /G S = 1.2 in Figs. 
10 (a) and (b), respectively. Fig. 10 (a) indicates that, for (G D /G S , ^,T) = (1.2,330,0), 
a metastable state appears at 5/j, = 0, while Fig. 10 (b) shows that, for (Gd/Gs, H,T) = 
(1.2,340,0), a stable state appears at intermediate 5/i. 




§H (MeV) Sn (MeV) 

Fig. 10: f2(T=0) against S/i for Gd/Gs = 1.2. The panels (a) and (b) correspond to \i = 330 and 
340, respectively. 



Numerical analyses show that, in the region where A = 0, Q decreases as 5/j, increases, 
so that is of the FM-type. Since A = at high temperature (cf. Figs. 3 and 4), Q is of 
the FM-type at high T. For the normal phase to appear, A ^ is necessary. 
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V. SUMMARY AND DISCUSSION 

In this paper, generalizing the analysis in |8j , we have dealt with unpolarized and polarized 
two-flavor quark matters at zero and finite temperatures. 

We have used the extended Nambu-Jona-Lassinio model and employed the mean-field 
approximation. The model contains two coupling constants, the quark-antiquark coupling 
constant G$ and the diquark coupling constant Gd- Although G$ and Go have the relation 
Go = (3/4)Gs, we have treated, as in [7|, Gd/Gs as a free parameter. 

The unpolarized case at T = are fully analyzed in [7|. In addition to it, we have 
analyzed the unpolarized case at finite temperatures through computing the thermodynamic 
potential. The result is summarized as the phase diagrams in the (p, T)-plane in Fig. 4. 
The structure of the diagrams are qualitatively the same for all values of Gd/Gs considered 
in this paper. 

We have dealt with two types of polarized two-flavor quark matters, i.e., the "magnetic- 
moment" -polarized and the "spin" -polarized quark matters. For describing such systems, 
the extended Nambu-Jona-Lassionio model is generalized to treat the chiral condensate, 
diquark condensate, and the degree of polarization on an equal footing. We have shown 
that, within the mean-field approximation, the form of the thermodynamic potential Q is 
the same for both the "magnetic- moment" -polarized and the "spin" -polarized cases. 

We have found that, as in the cases of other quantities, the S/j, dependence of Q heavily 
depends on the value of Gd/Gs- For small Gd/Gs < 1.15, VL is of the "ferromagnetic" - 
type, while, for large Gd/Gs > 1.15, there appears a window in the (/i, T)-plane, where 
Q is of the "normal" -type. We have found interesting behaviors of Q in the transition 
regions between the "ferromagnetic" -type and the "normal" -type regions. We have seen 
that, at high temerature and/or density, the "ferromagnetic" phase is energetically favored 
over the "normal" phase, which is in conflict with intuition. This fact may suggests that 
the extended Nambu-Jona-Lasinio model supplemented with the mean-field approximation 
is not applicable to the polarized quark matter in such regions. 

The analysis presented in this paper is the "first-stage analysis" . Closer investigation is 
to be done in various directions. Among others, are 1) inclusion of the residual interactions, 
and 2) incorporation of electron, positron, and gluons into the NJL model. Also interesting 
is to apply the method in this paper to quark matters in different phases, such as the gapless 
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charge-neutral 2SC phase, the color-flavor-loced phase, the non-uniformly condensed phases 
[lo| . and so on. 



APPENDIX A: THE ROLE OF THE PROJECTORS V T S (r = ±, s = ±) 



We start with briefly reviewing the spin content of q(x). We employ the box normalization 
by confining the system within a cube and introduce the periodic boundary condition to make 
the single-particle plane- wave basis. q(x) in the interaction picture is expanded, in standard 
notations, as 



q{x) 



E E \l ^§ MpK(p>~^ + &J03)i 



ip-x 



p s=± 



where V s (p)u s i(p) = 5 ss >u s (p) and V s {p)v s i{p) = 5 ss iv s (p). The third component of the spin 
operator reads S z = | / rf 3 xg^(a;)(T 3 g(a;). Expectation values of S* 2 in one quark- and one 
antiquark- states are computed, in respective order, as (cf. Q]) 

(0\a ± (p)S z ai(p)\0) 



2 
± 



-Tr P 



»(-)/ 



-p)a' 



(0\b ± (p)S z al(p)\0) 



m 3 
n 

? J 1 

2E r 



(A.l) 



u f (p)Pi(p)0- 3 :p ± (pMp) 



m 3 



"2£, 



(A.2) 



The minus sign on the right-hand side of the first line in Eq. ( I A. 21) comes from the fact 
that the quark obeys the Fermi-Dirac statistics. Then, we call the state with V+(p) = 1 
(V+(p) = —1) the "spin up" ("spin down") state. 

To see the role of P^, we first apply P« (iV) to q(x), 

P«(zV) ? (x) 



EE 




a s ,(p)V s (p)A4-p)u s/ (p)e- tp - x 



+bl(p)V s (-p)A-(p)vMe ip ' x ■ 
Noticing that A^(—p)u s /(p) = u s >{p) and A-(p)v s >(p) = 0, we have 



pi -) w) = E 




a s (p)u s (p)e~ 



ip-x 
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Thus P^(iV) projects out onto the "spin-up" ("spin-down") positive energy-state. Simi- 
larly, using A + (— p)u s '(p) = and K + {p)v s i{p) = v s >(p), we have 



Pi +) (^V) 9 (x) = £ \I^MvMe ip 




and then P^(iV) projects out onto the "spin-up" ("spin-down") negative energy-state. 



APPENDIX B: CONSERVATION OF Q + AND Q_ UP TO 0{G S ) AND 0(G D ) 

In this Appendix, we show that both the "magnetic- moment (MM) "-polarized net quark- 
number charges (Q + and QJ) and the "spin" -polarized ones (Q' + and Q'_) are conserved up 
to and including 0(G S ) and 0(G D ). 

The Hamiltonian reads 

H(x ) = J d 3 xH(x) =H + H S + Hps + H D , (B.l) 
H = q(-tf • V + m )q-G s [{qqf + (qhsrq) 2 ] , 

-G D [{ifee b l5 q) {iqee b l5 q°)] . (B.2) 

Here H , Hs, Hps, and Hp are, in respective order, the free-, the scalar quark- ant iquark in- 
teraction, the pseudoscalar quark- ant iquark interaction, and the diquark interaction Hamil- 
tonians. Since Q± = + and Q' ± = + <5±~\ it is sufficient to show that 
'H(x ),Q i :\x )'\ = 0(G s ,G 2 D ,G s G D ). 

The quark propagator G in the imaginary-time formalism is Gn | _ Q with m Q for m, 
where Gn is as in Appendix C: 

y- loP ( s ^ ("MM" -polarized case), 

(B.3) 



— 7oPs (p) — ( "spin" -polarized case) , 



where p = ipoE- 
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Derivation of Hq, Q 



We write Hn as 



H = J d 3 x q (— «7 • V + m ) q 

= Jd 3 xq^Ej2 [Pi _) («V) -Pi +) (zV)] q, 



where E = y — V 2 + m 2 ,. Using Eq. (12. 6p . we can easily show that the equal-time commu- 
tator between and H vanishes: 



[q { ;\h ] =o. 



Derivation of 



We start with computing 



[^s(ps), Qi r) ] = g s / rf 3 x [(g [pW(i v)r - rpW(iv)l g) (grg) 



+ (qTq) (q [p«(i V)r-rP«(iV) 



J g 

where T = 1 for and T = ry 5 r for -f/ps. 

Let us compute the ensemble average of Eq. (1B.4j) up to 0(Gs) using Eq. (1B.3[) : 

P0,90 

xTr {TPP(-p) - P< T) (g)r} G(p)r] , 



(B.4) 



(B.5) 



where V is the volume of the system. Appropriate regurarization is understood to be 
introduced. From Eq. (1B.3|) . we have 



(B.6) 



Using this relation in Eq. (1B.5j) . one can readily see that the right-hand side vanishes, so 



that 



H 



S(PS), Qs 



(r) 



0(G S , G D , G S G D ) 
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Derivation of H D , Q { I ] = 0(G 2 s ,G 2 D ,G D G s ) 



We start with computing 



[H D , QP] = 2G D J d 3 x[(q c ee b l5 P^(iV)q) (qee b l5 q c ) 
-(fee b l5 q) (gP«(i V)ee b l5 q c " 



Computation of the ensemble average of Eq. (1B.7j) up to 0(Gd) yields 

d 3 p f d 3 q 



xTr [G T (g)C 75 {P«(-£)G(d - G(p)pW(p)} 75 C] 



(B.7) 



(B.8) 



Using the relation (1B.6I) in Eq. (IB. 81) . one can readily see that the right-hand side vanishes, 



so that 



Hd, Q 



0(Gg, G 2 D , G S G D ) 



APPENDIX C: PROPAGATORS 

In this Appendix, for completeness, we write down the form of propagator G with G^ 1 
as in Eq. (13.101) . The propagator G for the blue quark is obtained from G by taking the 
limit A — > and deleting 1- 



Magnetic-moment-polarized quark matter 

Straightforward manipulation yields, with obvious notation, 

Po - E p =F H± s 



Cll(22) - ^2 



G 



21(12) 



;7o 



(p + E p ±fi TS )(p -E p TH±s)-\&\ 2 3 



Pi +, (P) 



(p -E p ±Li± s )(po + EpTfJ<Ts)-\A\ 2 
A± 



E 



(po + £ P ± p TS )(po ~ EpT M± s ) - |A| 2 
A± 



7oP ( . _) i 
P«(p) 



1/1 



(po - -Bp ± p± s ) (p + E p =f p TS ) - | A 



2 * 



where p = ipoE- 
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Spin-polarized quark matter 



The propagator G is deduced as 



G 



11(22) 



E 



Po - E p =f H- 



(p + E p ± fi B )(p -E p t H-,) ~ I A| 2 



7o 



Pl +) ( 



p + E p =f ji- s 



r7o 



(po-E p ±^)(p + EpTV-s)-W J s {P) _ 
^21(12) = E ~ ■ 5 i .. w_ A± r, — \ , A ., Pi +) (P) 



lfl 



(po + £ P ± A*s)(Po -^T li-,) ~ |A| 2 



(po -E p ± /U s )(p + ^ P =F A*- s ) - |A| J 



P(-)(P) 



APPENDIX D: DERIVATION OF EQ. (3.13) 



In this Appendix, we compute Det in Eq. (I37T2D with Eq. f[3~T0l) . 



Det = Det 



A + (Go)" 1 1/1 



Using the identities, 

Det ( A B \ = (-) n Det (CB - CAC^D) = (-)"Det (BC - BDB^A) , 

where A, B, C, and .D are (n x n) matrices, we obtain 

[Det {-(3G- 1 )} 2 = Det j/? 4 J] [a t A^ - A T (G^ 1 (A T ) _1 (G^T 1 ] J 

= Det j/? 4 J] ["|A| 2 - 75 (Gjr 1 75 (G^)" 1 



(D.l) 



where use has been made of Eq. (13. 7p . 

Substituting Eq. ([32]) into Eq. flD~2l) . and using Eqs. (jXT) . fl2T8j) . and (JZED we obtain, 
after some algebras, 

1 2 



[Det (-/3G" 1 )]' = Det [/3 4 {iq - 2p (p 2 - E 2 p - - |A| 2 ) 
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where po = ipoE (cf. Eq. (13. 9p ) and 

F 1 = (pi - El - - |A| 2 ) 2 + p 2 (//+ - /i_) 2 - E 2 p (//+ + /i_) 2 . 

Using Eq. (12. 5p with Eqs. (12. 3p and (12.41) . and Eq. (12. 2p . we obtain, after some algebras, 



[Det (-(3G- 1 )] 2 = Det£ 4 



-Fll -^12 



-^21 ^22 



Here 



11 — ^22 



F!l 2 + — (a-n)F 2 



2 J -/- L c » 



12 



F 21 = - — [(pxn)-fF 2 ] I/I" 



where n = n(p), 1 2 is the (2 x 2) unit matrix in a spin space, and 

F 2 = 2p - //_) (p 2 - £ 2 - - |A| 2 ) . 



Using Eq. (ID.ip again, we obtain, after some algebras, 

[Det (-(3G- 1 )} 2 = Det {/3 8 (F| - F 2 ) 1,1^} , 

where use has been made of \p x n| = p 2 n 2 — (p-n) 2 and Eq. (12.21) . Further manipulation 
yields (cf. Eq. ( Eg} ) 

In Z r>9 = 1 In Det (-(3G^) = N f £ £ In f J] J] [p 2 - (fig) : 



with 



E v + a^±±^) +|A| 2 



n p \P=± <7=± 

1 1/2 



+ P 



(D.3) 



Using the relation 



J> [/3 2 (pi - £ 2 )] =p[E + 2Tln(l + e^ E )] 



and making the replacement being valid in the large- V limit, 

d 3 p 



(2tt)' 
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we finally obtain for the contribution from the red- and green-quarks to the thermodynamic 
potential, 

In Z T 



n 



r, g 



-T- 



■"r,g 



V 



-Ni 



/t^E [« + 2Tln(l + e-^ 

J V I p,o=± 



(D.4) 
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